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We investigate the cosmological effects of neutrino lumps in Growing Neutrino Quintessence. The 
strongly non-linear effects are resolved by means of numerical N-body simulations which include 
relativistic particles, non-linear scalar field equations and backreaction effects. For the investigated 
models with a constant coupling between the scalar field and the neutrinos the backreaction effects 
are so strong that a realistic cosmology is hard to realize. This points towards the necessity of a field 
dependent coupling in Growing Neutrino Quintessence. In this case realistic models of dynamical 
Dark Energy exist which are testable by the observation or non-observation of large neutrino lumps. 


I. INTRODUCTION 

The origin of the observed accelerated expansion of the 
universe is still unknown mm- It is usually accounted 
for by a Dark Energy (DE) component. The simplest 
possibility consistent with observations is a cosmological 
constant A, but a lot of alternatives have been proposed 
[3]. Prime candidates are dynamical Dark Energy mod¬ 
els mediated by a scalar held or modihed gravity - the 
latter being often equivalent to the former [4] . For many 
alternatives the cosmological constant problem 0 IS] of 
explaining the small value of A persists, however. Also 
the explanation of why DE becomes important in the 
present cosmological epoch is often not more convincing 
than for a cosmological constant. 

Growing Neutrino Quintessence (GNQ) [71 |5] offers 
here some advantages. As a quintessence model 13 [10] 
the late time acceleration is driven by a scalar held ip (the 
cosmon), employing a mechanism similar to Inhation. It 
is possible to unify the late and early time acceleration 
into a single picture [IIHI3 so that the same held is re¬ 
sponsible for DE and Inhation. As an overall description 
within quantum gravity crossover cosmology |14j GNQ 
also addresses the cosmological constant problem. 

GNQ is able to explain the smallness of the DE com¬ 
ponent, since the dynamical DE density decays during 
the cosmic history, just as the other energy densities in 
the universe. The DE density being small is then just a 
matter of time - it is small because the universe is old. In 
contrast to simpler Quintessence models GNQ solves the 
Why-Now-Problem. No hne tuning of the self interaction 
potential is needed for this purpose. A coupling between 
the cosmon and the neutrinos provides a mechanism for 
stopping the evolution of the cosmon held as soon as the 
neutrinos become non-relativistic. The phenomenology 
of a very slowly evolving scalar held resembles a cosmo¬ 
logical constant. The transition from relativistic to non 
relativistic neutrinos acts as a trigger for the DE domi¬ 
nation. For neutrino masses allowed by observations this 
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transition happens in the “recent” past, explaining why 
DE has become important now. 

Despite a background evolution similar to the AGDM 
model for redshift z < 5, GNQ has a phenomenology 
which is distinct from other models. It predicts a time 
varying neutrino mass and the formation of neutrino 
lumps, which might be detectable through there grav¬ 
itational potentials m- The formation of lumps is a 
consequence of the large coupling between neutrinos and 
the cosmon, which is required for the stopping mecha¬ 
nism. The resulting additional attraction between neu¬ 
trinos is about 10^ times stronger than the gravitational 
attraction. It can have a natural explanation in a particle 
physics framework |S]. 

While the strong coupling on the one hand offers with 
the lumps a clear and distinct way of testing the model. 
On the other hand, it renders the model technically dif¬ 
ficult to study. In GNQ perturbations in the neutrino 
density become non-linear already at z « 1 — 2 on very 
large scales [15]. This has lead to the development of a 
comprehensive N-body simulation [mini to follow the 
formation of the neutrino lumps. The simulation is dif¬ 
ferent from the usual GDM-only simulations: In order 
to include backreaction effects, induced by the highly 
non-linear nature of the lumps |18j . the background is 
solved simultaneously with the perturbations. Addition¬ 
ally, neutrinos becoming relativistic during the formation 
of lumps is captured by the simulation. A similar frame¬ 
work for relativistic N-body simulation with focus on the 
metric perturbations was explored recently in [T 3 . With 
our simulation it was possible to draw a consistent pic¬ 
ture of neutrino structures within GNQ. Eor stable lumps 
the main characteristic features can be understood within 
an approximation in terms of a non-relativistic fluid of 
neutrino lumps [20] ■ 

In this work we investigate if GNQ can provide a real¬ 
istic expansion history. Therefore we study the equation 
of state and the energy density of the cosmon for different 
model parameters. We aim at finding model parameters 
for which the backreaction effect remains compatible with 
an accelerated expansion with Qde ~ 0.7. At the same 
time the accelerated expansion of the universe must start 
early enough to be consistent with observations. 
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A time dependent neutrino mass related to a scalar 
Dark Energy field concerns a wider setting than GNQ. 
Mass varying neutrino scenarios (MaVaNs) have been 
studied earlier in m and share common features with 
GNQ as the instability of neutrino perturbations [221l24j . 


This work is organized as follows. We start with a brief 
review of GNQ in section [IT] In section m we discuss 
the formation of lumps and their backreaction on the 
cosmological expansion. In section IV we describe our 
simulation, which we use to perform a parameter scan. 
Results are presented in section |Vj Finally, we conclude 
in section ED 


II. GROWING NEUTRINO QUINTESSENCE 
A. Basic Concepts 

In this section we briefly describe GNQ. The ingre¬ 
dients of GNQ are a scalar field (p (the cosmon) and 
neutrinos. The neutrino mass depends on the value of 
ip, thereby coupling the cosmon and the neutrinos. The 
cosmon itself is described by the standard Lagrangian 
of a scalar field which takes, using the metric signa¬ 
ture (—,-f,-f,-f) and setting the reduced Planck Mass 
to unity, SttG = 1, the form: 

+ V{p). (I) 

We choose an exponential potential V{p) oc As 

long as the neutrino mass can be neglected the expo¬ 
nential potential leads to scaling solutions of the cosmon 
field. The background energy density of the cosmon be¬ 
comes independent of the initial conditions and mimics 
matter (radiation) during matter (radiation) domination 
[5] , where the energy density of the cosmon is a constant 
fraction of the total energy density Here w 

is the equation of state of the dominating species. Gon- 
straints on early dark energy (EDE) require a > 10 [25l - 
I28j , where we use a conservative bound in view of possible 
unexplored parameter degeneracies. 

The dependence of the neutrino mass on the cosmon 
is given by: 


/3 = 


dp 


( 2 ) 


In general the coupling /3 can be (/^-dependent. We estab¬ 
lish in this note that the size of the backreaction effect 
depends crucially on the presence or absence of a varia¬ 
tion of P{p)- An investigation of a particle physics mo¬ 
tivated variation of [5] in reference has revealed a 
small backreaction effect and an overall cosmology con¬ 
sistent with present observations. For a constant /3 large 
backreaction effects have been observed [16]. We adress 
here the question if the model remains compatible with 
observations in this case as well. 


A constant coupling implies for the neutrino mass: 


m,,{p)=m^e (3) 

where an additive constant in p is fixed such that V{p = 
0) = 2.915 • 10“^ eV. The (^-dependent neutrino mass 
allows for energy transfer between neutrinos and the cos¬ 
mon, which is proportional to the trace of neutrino energy 
momentum tensor: 

= +PT^^^p 

= -I3T^.)P. (4) 

The trace of the energy momentum tensor = 

, = —p^ -f 3Py vanishes for ultra-relativistic neutri- 

nos. The coupling between neutrinos and the cosmon is 
therefore ineffective for relativistic neutrinos. The neu¬ 
trino energy-momentum tensor also sources the Klein- 
Gordon equation which governs the evolution of the cos¬ 
mon: 


V^V^p-V'ip) = PT^,y (5) 

We will describe neutrinos and dark matter by an N-body 
simulation. The trajectories of classical neutrinos obey a 
modification of the geodesic equation m- 

ci?7^ 

— + = Pd^^p + Pu-'d^pu^^, ( 6 ) 

dr 

where denotes the four-velocity and r the proper time. 
The left hand side is the usual gravitational motion, with 
the Ghristoffel symbols T^,^ determined by the metric. 
Throughout this work we use the Newtonian gauge for 
the metric: 

ds^ = - (1 + 2di) df + {1 - 2^) . (7) 

We will work to first order in the gravitational potentials 
$ and T and neglect there time derivatives. 

The right hand side of equation § describes an addi¬ 
tional force due to the coupling to the cosmon. It consists 
of two parts. First, a velocity dependent part f5u''dvpu^ 
compensates changes in the mass for neutrinos moving in 
a varying cosmon field so that momentum is conserved. 
A neutrino moving into a region with smaller (larger) 
values of p will lose (gain) mass. To compensate the loss 
(gain) of momentum it will be accelerated (decelerated). 
Second, a velocity independent fifth force I3d^p. In the 
non-relativistic limit it acts as an attractive force about 
2/3^ times stronger than gravity |2D|- 


B. Homogeneous Evolution 

Let us now turn to the homogeneous limit and discuss 
how GNQ in its simplest form can lead to an acceler¬ 
ated expansion of the universe. At early times when the 
neutrinos are relativistic the evolution of the cosmon is 




3 


p 


Mp^ Hq^ 



FIG. 1. Energy density of the cosmon-neutrino fluid, for 
different parameters a and j3. We compare to the CDM- 
density and the density of a cosmological constant Ha- The 
parameters where chosen to match Q,a today. The stopping 
occurs earlier for larger a, with a smaller amount of early 
dark energy. 


determined by the potential. Therefore the cosmon will 
evolve towards its scaling solution with the DE density 
decreasing with a~^ during matter domination. In view 
of the growing mass the neutrinos become non-relativistic 
rather late. The interaction becomes important once 
aV(ip) « « —fipu- It acts as an effective potential 

barrier stopping the time evolution of the energy den¬ 
sity of the cosmon-neutrino fluid. The constant energy 
density then mimics a cosmological constant. Since the 
energy density of the neutrinos is small compared to the 
cosmon energy density, the coupling must be rather large. 

Most of the cosmological parameters as Ude = 
and mu are approximately independent of the individual 
values of a and /?, they only depend on there ratio. De¬ 
manding a dark energy density of IIde ~ 0.7 enforces 
— £ ~ 5 [7] for a present neutrino mass = 0(1 eV), 
where smaller neutrino masses require large — We 
note that the usual cosmological bounds on the neutrino 
mass from CMB and Large Scale Structure observations 
[301 El] do not apply here, since neutrino masses have 
been substantially smaller in the past. In the homoge¬ 
neous limit the neutrino mass is mainly constrained by 
earth based experiments. Also the scale factor at which 
the neutrinos stop the cosmon evolution has only a mod¬ 
erate dependence on the individual values of a and j5. 
The energy density fraction of the cosmon before stop¬ 
ping is given by cx a~‘^ and hence becomes smaller 
for larger a. The time at which the interaction with 
neutrinos compensates the self interaction of the cosmon 
becomes earlier for larger a. The onset of dark energy is 
therefore earlier for larger values of a and 13. 

As we will discuss later, strong backreaction effects 
will alter this simple picture. We will see in section |III| 
that backreaction effects always counteract the stopping 



FIG. 2. The scale factor onl at which the dimensionless linear 
neutrino power spectrum becomes unity, A(fc, anl) = 1 , as a 
function of scale, for the parameters a = 10 and j3 = —52. 
Already at a ~ 0.40 scales around k ~ 0.02 are non-linear 
demonstrating the failure of standard perturbative methods 
compare to the same figure in reference m- 

mechanism and the cosmon will evolve again, so that it 
is not guaranteed that values for a and f3 which describe 
a realistic cosmology in the homogeneous limit will also 
describe a close-to-realistic cosmology including backre¬ 
action. 

Since backreaction effects can only be important af¬ 
ter the neutrinos became non-relativistic the homoge¬ 
neous description remains valid at early times. Large 
values for a are preferred by bounds on early dark en¬ 
ergy. For large a the stopping mechanism acts earlier, 
hence also the backreaction becomes important earlier. 
From these qualitative considerations we already find 
some tension between reducing the backreaction effects, 
which spoil the stopping of the cosmon evolution, and 
satisfying bounds on EDE. 


III. BACKREACTION AND EFFECTIVE 
EQUATION OF STATE 

A. Neutrino Lumps 

In GNQ it is important to understand structure for¬ 
mation, not only in view of using Large Scale Structure 
observation as a probe for our cosmological models, es¬ 
pecially to test DE models or “measure” the neutrino 
mass. It is crucial to understand the formation and evo¬ 
lution of neutrino lumps, before being able to judge about 
the viability of GNQ as a DE model. In this section we 
shortly review the progress towards an understanding of 
the neutrino lumps, for details we refer to previous work 
[m US m Eoi m EH EiEi ■ Our main focus lies on the 
strong backreaction effects from non-linear perturbations 
in the neutrino-cosmon fluid. 

The large non-linearities have there origin in the large 
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coupling P = O(IO^). Therefore the additional force be¬ 
tween neutrinos will be about 10^ —10^ times larger than 
the gravitational interaction between neutrinos and be¬ 
tween neutrinos and CDM. In turn the neutrino pertur¬ 
bations grow very quickly as soon as neutrinos become 
non-relativistic. This implies that the fluctuations in the 
neutrino energy density become non-linear even at large 
scales. The scale factor unl at which this happens for a 
neutrino perturbation of a given wavelength k~^ can be 
estimated by the value of a at which the linear dimension¬ 
less neutrino Power Spectrum A^{k) = fc^Pi,(fc)/(27r^) 
becomes order unity. Looking at figure we see that for 
the particular choice of parameters a = 10 and /3 = —52 
already at a ^ 0.4 scales around /cnl.i/ 0.01 hMpc“^ 
become non-linear, while today scales around 
0.002 h Mpc“^ are non-linear. The exact value of the 
non-linear scale of neutrino-cosmon perturbations de¬ 
pends on the chosen parameters, but it is a generic find¬ 
ing that is smaller than the corresponding wave 

vector for CDM perturbations, A:nl,c,o ~ 0-1 hMpc“^. 
These can be traced back to instabilities in the neutrino 
perturbations already present at linear order. These in¬ 
stabilities are stabilized non-perturbatively by the forma¬ 
tion of neutrino lumps. 


B. Backreaction 

Usually backreaction in cosmology is assumed to be 
negligible. In the last years several quantitative esti¬ 
mates [551 - I57| came to the conclusion that backreaction 
is indeed small in the ACDM-model. In contrast, backre¬ 
action effects are crucial in GNQ. We demonstrate this in 
figure where we compare the numerical results for the 
clumping neutrinos with the pure background evolution 
for which the effects of non-linear neutrino perturbations 
are neglected. We choose the parameters a = 10 and 
/3 = — 52 that have often been employed in the litera¬ 
ture. 

We find two types of backreaction effects. First, the 
Friedmann equation involves the volume averaged energy 
density, which we will dehne below. Second, the average 
value of the cosmon (f can not be obtained by solving the 
homogeneous equation of motion. The Klein-Gordon- 
Equation needs to be modified to include backreaction 
effects from the neutrino lumps. The reason is that the 
typical velocities and masses of the neutrinos do not co¬ 
incide with there counterparts of the homogeneous calcu¬ 
lation. While the first effect mainly affects the expansion 
history of the universe, the second effect is also important 
for an understanding of the lump dynamics. 

Let us first discuss the second effect. Due to the strong 
interaction most of neutrinos are bound in the lumps. 
Inside gravitational bound objects the gravitational po¬ 
tential has a well. Similar, inside neutrino lumps the 
local field value is smaller than its average by an amount 
of Sip. The mass of a neutrino inside a lump is there¬ 
fore smaller than the mass of a free-streaming neutrino 


^DE 




FIG. 3. Dark energy density fraction Ude (top) and equation 
of state w (bottom) as a fnnction of the scale factor, for a = 10 
and /3 = —52, with and without backreaction. 

m{(p + Sip) < m{<p) . As a consequence most of the neu¬ 
trinos have a mass substantially smaller than the mass 
estimated from the homogenous calculation. Due to the 
velocity dependent force the loss of mass during the for¬ 
mation of lumps is accompanied by an acceleration to rel¬ 
ativistic velocities. These two effects lead to a mismatch 
between the energy momentum tensor of neutrinos from 
the homogeneous calculation and its average value, as 
soon as the formation of lumps has started. 

We account for the backreaction effects by using the 
volume averaged energy momentum tensor. The Klein- 
Gordon equation for the average held is given approxi¬ 
mately by: 

0 + 3H0 + aVi<p) = -PT^,), ( 8 ) 

where the volume average is dehned as 

3$)r(,). 

(9) 

The determinant of the spatial 3-metric up to hrst 
order in metric perturbations is given by \/g^ ~ 
a? (1 — 3$). The integration is to be understood over 
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the whole simulation box. The volume is given by 
V Ki J (fix (1 — 3$). Taking backreaction effects con¬ 
sistently into account and evolving the volume averaged 
field (p additional modifications arise in the equation. 
However, we will neglect these terms for the qualitative 
discussion of backreaction in this section and postpone a 
more detailed discussing to section 

The right hand side of equation (^can be written as: 

= j3 {-pi, + 3P„) = -j3pu (1 - 3^,,) < -jipu, 

( 10 ) 

where the energy density and pressure are understood as 
volume averages. We use them to define the equation 
of state w,j. The neutrino pressure is positive > 0) 
such that pressure effects lower the effective potential 
barrier which stops the cosmon evolution. As a conse¬ 
quence, the time at which the cosmon evolution stops is 
postponed towards the future. If the evolution has al¬ 
ready stopped the effective reduction of the barrier can 
have the effect that the cosmon will evolve again. The 
weaker interaction between the neutrinos and the cosmon 
after the formation of lumps, can also be interpreted as 
a lower effective coupling /3;, which gets renormalized by 
integrating out short wavelength modes [20]. In a quali¬ 
tative sense f5i can be interpreted as the effective coupling 
between a fluid of neutrino lumps and the homogenous 
cosmon field. The smaller value of j3i as compared to (3 
is the dominant backreaction effect in our model. 

We next turn to the backreaction effect for the evo¬ 
lution of the background metric. One needs to replace 
the background density of neutrinos and the cosmon by 
their volume average, such that the Friedmann equation 
becomes: 

= Pc^DM + pv + pip. (II) 

In the presence of lumps pu has contributions from the 
neutrino velocities, and p^p involves additional gradient 
contributions. The observable DE component is the com¬ 
bined neutrino-cosmon fluid pde- The neutrinos are typ¬ 
ically subdominant but still contribute a significant frac¬ 
tion ^ 0.1. With an equation of state Wi, ~ 0.1 
the neutrinos lift the dark energy equation of state away 
from iti « — 1 to some higher value. 

The volume average of the cosmon energy density is 
given by: 

P<P = {djip)5^3 -h !/(</?), (12) 

where we only keep metric perturbations up to first or¬ 
der, neglected their time derivatives and use that the 
volume average of the gravitational potentials vanishes 
= T = 0. Also assuming that time derivatives of the 
cosmon perturbation 5ip are small allows us to approxi¬ 
mate « <p'^. Using the quasi static approximation is 
justified although the individual neutrino velocities are 
large. For the quasi static approximation to hold it is 


sufficient that the energy-momentum tensor for all neu¬ 
trinos does not evolve fast, so that there are no fast vary¬ 
ing sources for the cosmon. A non-zero Sp results in 
a positive contribution to the pressure, making it even 
harder to achieve an almost constant energy density for 
the cosmon-neutrino fluid. 

Without the gradient term one has the usual competi¬ 
tion between potential and kinetic energy. The potential 
energy should be dominant in order to have an accel¬ 
erated expansion. The averaged potential energy V{ip) 
differs from the potential energy V{0) of the averaged 
field p only by a few percent, such that no major back- 
reaction effect arises from this source. In contrast, the 
gradient term can be almost as large as the the potential 
energy. From the expression for the pressure 

Pv ~ ~ 

we see that a gradient term dominated equation of state 
would be Wjy = — ^. We emphasize that all backreaction 
effects individually lead to an evolving energy density of 
neutrino-cosmon fluid and typically push w away from 
- 1 . 

For models with constant /3 the lumps have the ten¬ 
dency to stabilize and to remain present once formed. 
The neutrino-cosmon fluid can be understood as an ef¬ 
fective fluid of nearly virialized neutrino lumps with pa¬ 
rameters differing from the microscopic ones |2()j . The 
observable DE is then the sum of a neutrino lump fluid 
and a homogenous background field. For virialized lumps 
the pressure between relativistic neutrinos and cosmon 
gradients is expected to cancel m- Therefore the equa¬ 
tion of state of the lump fluid is close to zero, similar 
to the fluid of non-relativistic neutrinos. The backreac¬ 
tion effect that remains even in this limit is the reduced 
effective coupling /3; between neutrino lumps and the cos¬ 
mon background field. Due to the not completely virial¬ 
ized lumps the pressure contribution from the neutrinos 
and the cosmon gradients do not cancel exactly, adding 
a small but relevant additional backreaction effect. This 
is different to gravitationally bound objects, for which 
a non-renormalization theorem states that small virial¬ 
ized objects decouple completely from the background 
evolution and there is no backreaction effect from small 
virialized objects at all [55] . 


IV. N-BODY SIMULATION 

The highly non-linear nature of the neutrino lumps 
makes there description non-amenable to standard per¬ 
turbative techniques. Instead we use a N-body simula¬ 
tion specially designed for GNQ. The N-Body simulation 
solves the background and the inhomogeneities simulta¬ 
neously and therefore allows us to study the backreaction 
effect of lumps on the homogeneous background evolu¬ 
tion. Concept and many details of the simulation were 
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already described in [Mini, we focus here on the equa¬ 
tion of motion for the average cosmon field (p and its 
perturbation bp. 

In our simulation we follow the usual motion of non- 
relativistic CDM particles and there clustering due to 
gravity. In contrast to the standard picture of structure 
formation the two gravitational potentials differ, $ 7^ 
because of the anisotropic stress from the neutrinos. 
This is accounted for by solving the Poisson equation for 
d) — 'h, which yields $ ones the Newtonian potential 'h 
is known. The Poisson equation for T is sourced by the 
energy density of CDM, neutrinos and to a small part by 
the one of the cosmon perturbations. 

The neutrinos are evolved using equation ([^. The 
cosmon evolution is governed by the Klein-Gordon equa¬ 
tion 1^. We split the cosmon into the volume average 
P = y f d?x \/g^(p and a perturbation 5ip = p — (p. Ne¬ 
glecting time derivatives of the gravitational potentials, 
time derivatives commute with the process of averaging 
(p K p. The averaged equation (H is: 

(p + ‘^H(p + a{l + 2'^)V{(f) 

= -/3(l + 2^)T(,) + (14) 


where we expanded up to first order in metric perturba¬ 
tions. Equation (14) is the full version of equation ([^ . As 
already discussed in section |III| the most important dif¬ 
ference as compared to a naive homogeneous calculation 
is the use of the actual average of the neutrino momen¬ 
tum tensor. Including the gravitational potential in the 
average gives only a minor correction. Also the averaged 
potential term agrees up to a few percent with the ho¬ 
mogeneous estimate. The gradient terms is roughly one 
order of magnitude smaller than the potential term and 
therefore only subdominant. Nevertheless, our numerical 
code includes all these effects. 

By subtracting equation (14) from the the Klein- 
Gordon equation © we find the equation for the per¬ 
turbation: 


Sip + SHSp — a didj5p{l + 2^) 

- $)) {dip) + {dip) 

-ba ((1 -f 2'i>)V{p) - (1-b 24')P((^)) 

= -/3((l + 2vI-)T(,)-(rT2^)^). (15) 


This equation is a non-linear wave equation, which is, 
due to the averaging, non-local in position space. To be 
able to solve this equation we need to make some ap¬ 
proximations. We employ a quasi static approximation 
for the cosmon perturbation for which we neglect the sec¬ 
ond order time derivative dp. Simply neglecting all time 
derivatives is not a consistent approximation. Doing so 
the resulting equation does not ensure that the perturba¬ 
tion has a vanishing mean 6p = 0. This can be seen by 
averaging equation (15). Taking into account the $ de¬ 
pendence in the volume element and only keeping terms 


to first order in the metric perturbations all terms except 
the time derivatives cancel: 

Sp + 3H6p = 0. (16) 

This relation ensures that the if the average vanishes ini¬ 
tially it will vanish at all times. This is still true if we 
neglect the second time derivative while keeping the first 
one. This approximation is consistent with the approxi¬ 
mation for kinetic term of the average energy density and 
pressure 


(^2 _ -b 

where we neglected the dp’^-teiui. So we neglected all 
terms which are second order in the time derivatives of 
the cosmon perturbations these terms are smaller than 
those with only one time derivative. 

If one instead neglects the second derivative with re¬ 
spect to conformal time the Hubble damping changes 
3H — 2H, we compared both possibilities and found 
only a small difference. We interpret this a sign that the 
quasi static approximation is justified. 

To solve the equation for 5p we use a Newton-Gaufi- 
Seidel (NGS) multigrid relaxation method, already ap¬ 
plied to the varying coupling model m and originally 
developed for modified gravity [38]. The quasi static 
approximation is crucial for applying the NGS method, 
which is not applicable to wave like equations, but can 
be applied to diffusion like equations [35]. The idea of 
the NGS solver is to rewrite the equation to be solved 
into a functional form: 


C[5p\ = DSp — F[Sp] = 0, 


(18) 


with some differential operator D and a non-linear func¬ 
tional F. The root of C[5p\ = 0 can be obtained by a 
newton-like iterative procedure: 




-1 


(19) 


the derivative is taken at each point individually, the 
coupling between different points, induced by the deriva¬ 
tives, is taken into account solely by the iterative pro¬ 
cedure. The derivative of the differential operator 
is defined by the discretisation rule used in the simula¬ 
tion. We define the gradient and the laplacian by relating 
a grid point to its neighbors in j-direction by a Taylor 
expansion: Sp{xi ± AxSij) = 5p(xi) ± djSp{xi)Ax + 
^dj6p(xi)Ax‘^ -b ..., with Ax the spacing between two 
grid points. The laplacian is then approximated by a 
seven-point stencil and the derivative is —6/Aa:^. The 
derivative of the gradient vanishes. 

In principle this method can be applied even in the 
presence of the non-local terms present in equation 
(15 ).In practice this not possible because calculating the 


non-local terms involves an integration over the full sim¬ 
ulation box in each iteration step. We account for these 
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terms iteratively. The difference between the values of 
the average terms of two time steps is small. So we use 
at a given time step the average terms of the proceed¬ 
ing time step as first approximation and apply the NGS 
solver a few times to correct for the difference. 


V. RESULTS AND DISCUSSION 

Using the N-Body simulation described in section [Tv| 
we perform a parameter scan and search for parame¬ 
ters describing a realistic universe with accelerated ex¬ 
pansion. For the details on the formation of lumps and 
there characteristics we refer to previous work [Milo]- 
We use a simulation box with a comoving volume of 
V = (600 h~^Mpc)^, which we divide into Nc = 128 
cells. The number of effective CDM particles Nq and 
neutrino particles N,^ is chosen to be equal to the num¬ 
ber of cells Nc = Nc = N^. The initial power spectrum 
has a spectral index of Us = 0.96 and an amplitude of 
As = 2.3 • 10“® at the pivot scale fcpivot = 0.05 Mpc“^. 
We start our simulation with the CDM particles only at 
flini.c = 0-02 and add the neutrinos at a later time, after 
they became non-relativistic. 

In view of the strong backreaction effects it is no longer 
clear that the stopping power of neutrinos for the time 
evolution of the cosmon is sufficient in order to account 
for a large present fraction of dark energy and an acceler¬ 
ation of the expansion similar to a cosmological constant. 
If so, the parameter range where this happens may be 
rather different from the one where the background evo¬ 
lution neglects the effect of neutrino structures. 

Our model has three parameters relevant for this in¬ 
vestigation, namely a, related to the amount of EDE, the 
neutrino cosmon coupling /3 and rrii, related to the size of 
the neutrino mass. We have performed a parameter scan 
in order to search for a parameter range consistent with 
observations. For this purpose we vary the parameters 
a and /3 individually while fixing the mass parameter to 
rrii = I eV. Figure shows that changing the mass pa¬ 
rameter by a factor of 10 effects the effective equation of 
state and the energy density by no more than 10%. 

A realistic DE model must certainly assume the bench¬ 
mark values for the present DE density Dde.o ~ 0.7 and 
the present equation of state Wq ~ — 1- In figure we 
show the values of Ude.o and wq for a grid in the param¬ 
eter space for a and /3. Sufficient acceleration typically 
requires rather small values a ;< 5. A band with an ac¬ 
ceptable fraction of present DE is typically found in the 
range 5 a 10, showing some tension already at this 
stage. 

The parameter range yielding an accelerated expansion 
{a < 5) is problematic also in view of the bounds on EDE, 
which require a > 10. In the parameter range where 
one finds wq < —0.9 some tension persists if one tries 
to get both the equation of state and the energy density 
compatible with observations. For a = 3 and a = 4 we 
indeed find wq < —0.9 but the energy density exceeds 


^DE 




FIG. 4. Energy density fraction of the cosmon-neutrino fluid 
Dde and effective equation of state w, for different mass pa¬ 
rameters mi, with a = 5 and /3 = —78. Even for mass param¬ 
eters different by a factor of 100 the equation of state varies 
at maximum about 10%, indicating that the value of plays 
only a minor role. 


with Ode ~ 0.75 the benchmark value of Ode,o ~ 0.7. 
On the other hand for a = 5 one has Ode ~ 0.7, but 
the equation of state is wg « —0.7. Although we could 
not find parameters for which wq and Ode,o rnatch the 
benchmark values precisely, our results are not too far 
from those values either. It might be that varying also 
the mass parameter rui could bring them into agreement 
with observations. 

The equation of state is not constant in time, it can 
even possess oscillating features, see figure It may 
happen that the present time coincides with a minimum 
(maximum) of w during an oscillation. In this case the 
cosmic evolution is actually better described by an aver¬ 
age value somewhat larger (smaller) than wq. The time 
evolution of the equation of state is shown in hgure 
for a range of parameters a and j3 in the region not too 
far from the benchmark values. One typically observes 
a first stop of the scalar field {w « —1). Due to back- 
reaction this is followed by a slow decrease of the dark 
energy, typically with —0.9 ^ w < —0.8. 

Only looking at the energy density and the equa- 
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FIG. 5. Present energy density SIdb.o and equation of state 
Wo of the cosmon-neutrino fluid. Realistic values (wo ~ — 1, 
flaE.o ~ 0.7 ) are found for small values of a. It is hard to get 
both values “correct” simultaneously, for sufficiently large a. 


tion of state today is not sufficient. In the parameter 
range acceptable for the benchmark the neutrinos be¬ 
come non-relativistic late. Consequently the cosmon evo¬ 
lution stops late. This is visible in figure where the 
first pronounced minimum in w precisely corresponds to 
the time when the increase of ip is first stopped and the 
oscillations set in. Supernova observations probe the ex¬ 
pansion history up to redshifts higher than z « I and 
prefer an almost constant Dark Energy |40] . We find that 
for close-to-realistic models the equation of state reaches 
values around —I only for scale factors a > 0.6, which is 
difficult to get into agreement with m — 1 from a ^ 0.5 


W 



--- a=4,yS=-52 - a=3,j3=-52 

FIG. 6. Equation of state as a function of the scale factor. The 
model parameters are chosen such that w and Ude.o are near 
the benchmark values. Values wo < —0.9 are only reached be¬ 
fore backreaction effects become important. Thus w ~ —0.99 
for 0 = 5 and /3 = —52 is not accompanied by large negative 
w at redshifts relevant for supernova observations. 


until today. 

Figure shows the generic evolution of the equation of 
state: It drops down after the neutrinos became non- 
relativistic followed by a few damped oscillations. In 
the homogeneous evolution these oscillations are damped 
away quickly and the equation of state assumes an al¬ 
most constant value rather close to w = — I. In fact the 
equation of state grows again due to the backreaction 
and typically reaches values w « —0.8. An equation of 
state of Wq < —0.9 is only reached before or shortly af¬ 
ter backreaction becomes important. This simply means 
that lumps had not enough time to grow large enough 
for being able to induce significant backreaction effects. 

From these results we conclude that GNQ with a con¬ 
stant coupling (3 is probably not a viable DE model. Re¬ 
alistic values for wq and Ude.o seem only possible if the 
cosmon evolution is stopped late, so that backreaction 
effects have no time to become important. Stopping the 
cosmon evolution late is in some tension with supernova 
data and involves a large amount of EDE, probably not 
consistent with observations. 


VI. CONCLUSION 

We have performed a numerical analysis of Growing 
Neutrino Quintessence with a constant cosmon-neutrino 
coupling /3. Due to strong backreaction effects from the 
formation of large neutrino lumps these models have dif¬ 
ficulties to be compatible with the observed properties of 
dark energy. 

A specific choice for the model parameters a, /3 and 
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rrii, which appears to be compatible with observations 
at the homogenous level, is typically no longer viable if 
backreaction is included. Our parameter scan reveals re¬ 
gions for which the backreaction effects are small enough 
to allow a slowly evolving cosmon and consequently an 
almost constant DE density. However, this is only possi¬ 
ble if the neutrino lumps form late so that backreaction 
effects are still small today. In this case an accelerated 
expansion is only possible for scale factors a > 0.6, in ten¬ 
sion with an almost constant equation of state for scale 
factors a < 0.5, as preferred by supernova data. Fur¬ 
thermore, the parameter region for which the equation 
of state is close to —1 and the DE density is not too far 
from 0.7, requires a < 5. This contradicts constraints on 
early dark energy for which a > 10 is necessary. We con¬ 
clude that growing neutrino quintessence with a constant 
coupling /3 is probably not a viable DE model. 

These results for a constant coupling should be con¬ 
trasted with models where (3 increases with Lp. For this 
second class of models the backreaction effect is found 
to be small since the neutrino lumps form and disrupt 
periodically m- At the present stage this second class 
of models seems compatible with observations. In cer¬ 
tain parameter ranges it may even be hard to detect a 
difference from the ACDM models and its variants. 

These two classes of models may be seen as particular 
points in a larger class of models where (3 is allowed to 
vary with (p. Having established points that are viable 
with only rather small deviations from ACDM, as well as 
other points where the deviations are so strong that the 
model is no longer acceptable, we can conclude by conti¬ 
nuity that in between there will be models which are still 
compatible with observations today, but also offer highly 
interesting prospects of finding deviations from ACDM. 
Finding large neutrino lumps, thereby observing the cos¬ 
mic neutrinos directly, would be a direct hint for GNQ. 
Even for models with small neutrino perturbations we ex¬ 
pect observable deviations from the ACDM model, due 
to the different evolution of the neutrino sector. First, 
the transition of relativistic to non-relativistic standard 
massive neutrinos is imprinted in the CMB fluctuations 
as well as in the matter distribution, with a specific scale 
dependence [sain]- The signal differs for constant or 


time-varying neutrino masses. Second, free-streaming 
standard massive neutrinos attenuate the growth of mat¬ 
ter perturbations on small scales and therefore add an 
additional scale dependent effect to the matter distri¬ 
bution. Observing these scale dependent effects as pre¬ 
dicted for standard neutrinos with a constant mass would 
be a strong argument for the ACDM model and against 
GNQ. 

The result for models with constant (3 presented in 
this note as well the results on the varying /3 model pre¬ 
sented in m suggest that only those models are viable in 
which the small scale non-linear neutrino perturbations 
have only a moderate effect on the large scale dynamics. 
Nevertheless, the neutrino lumps can have a observable 
effects on larger scales. One possibility to account for 
these effects is to construct an effective fluid for the long 
wavelength perturbations by averaging over small scales 
non-linearities as proposed in reference [36]. A similar 
route has already been taken in [20] to describe the large 
scale dynamics of virialized neutrino lumps in the con¬ 
stant (3 model by means of an effective lump fluid. These 
ideas where already successfully applied to the mildly 
non-linear regime of structure formation in the form of 
the Effective Field Theory of Large Scale Structure gT]- 
see also m- Adopting these ideas to GNQ we hope 
that it will become possible to study the dynamics of per¬ 
turbations in GNQ on large scales qualitatively. It might 
even become possible to study some effects of lumps on 
the CMB, without running time consuming simulations. 
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